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This  paper  presents  the  results  of  the  investigations  of  the  thermal  behavior  of  a  model  valve-regulated 
lead-acid  cell  at  overcharge  and  saturation  level  around  90%.  The  heat  capacity  and  the  heat  transfer  coef¬ 
ficient  are  experimentally  obtained.  The  quasi-steady  state  characteristics  of  the  cell  at  constant  charging 
current  are  measured.  The  experimental  data  show  that,  when  operating  in  recombination  regime,  the 
stationary  cell  temperature  is  a  linear  function  of  the  charging  current.  When  the  cell  is  gassing,  part  of 
the  heat  is  dissipated  with  the  released  gases,  as  a  result  of  which  the  cell  temperature  is  lower  than 
predicted  from  the  linear  temperature/current  dependence.  Introduced  dimensionless  slope  of  the  tem¬ 
perature/current  dependence  is  a  characteristic  of  the  VRLA  battery,  which  could  be  used  for  fast  and  easy 
estimation  of  the  temperature  increase  of  the  battery.  The  kinetic  parameters  of  the  oxygen  evolution 
reaction  (OER)  proceeding  at  the  positive  plates  of  VRLA  cells  are  experimentally  determined.  The  cur¬ 
rent/overpotential  dependence  of  the  OER  obeys  the  Tafel  relation.  The  temperature/current  dependence 
and  the  kinetics  of  the  OER  are  used  to  estimate  the  critical  charging  voltage  for  safe  operation  of  the 
battery. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  thermal  behavior  of  lead-acid  batteries,  and  especially  of 
their  recombinant  versions,  receives  considerable  attention  in  the 
context  of  its  contribution  to  the  heat  generation  in  power  systems, 
efficiency  of  operation,  reduction  in  service  life  and  susceptibility 
to  thermal  runaway.  The  thermodynamics  of  the  primary  reactions 
in  a  lead-acid  cell  is  described  in  [1].  The  thermal  “vicious  cycle” 
of  vented  batteries  with  low  heat  transfer,  caused  by  secondary 
reactions  of  water  electrolysis  during  overcharge,  is  analyzed  in  [2]. 
Thermal  runaway  is  a  result  of  the  positive  feedback  between  the 
heat  generation  due  to  Joule  losses  and  the  rate  of  Arrhenius  gas 
evolving  reactions  during  water  electrolysis. 

The  heat  generation  in  recombinant  VRLAB  on  overcharge  is  dis¬ 
cussed  in  [3],  where  in  addition  to  electrolysis,  the  internal  oxygen 
cycle  in  which  the  full  power  of  the  charging  unit  is  used  only 
for  heating,  is  considered.  The  VRLA  battery  on  overcharge  can 
be  generally  treated  as  parallel-connected  electrolyser  and  oxygen 
generator/recombinator.  The  thermal  stability  of  the  VRLA  battery 
is  determined  by  its  reduced  thermal  mass  and  thermal  conduc¬ 
tivity  of  the  interior  as  part  of  the  electrolyte  is  replaced  by  gas, 
both  decreasing  with  drying  [4,5].  Stability  analysis  of  VRLA  batter¬ 
ies  subjected  to  abusive  tests  in  high  temperature-high  overcharge 
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voltage  environment  is  performed  and  the  obtained  results  are  dis¬ 
cussed  in  refs.  [5-8].  The  influence  of  a  preceding  discharge/charge 
cycle  on  the  thermal  instability  is  demonstrated  in  [4,9,10,11  ].  Inves¬ 
tigations  of  the  thermal  phenomena  have  been  developed  in  several 
areas  as  thermal  mapping  [7]  to  visualize  the  temperature  dis¬ 
tribution  inside  the  cell;  thermal  video  imaging  and  calorimetric 
measurements  [12-14]  to  study  the  thermal  problems  arising  in 
EV  and  HEV  batteries;  several  electrochemical  studies  have  been 
performed  to  elucidate  the  influence  of  separator  structure  [15], 
saturation  level  [16,17]  and  electrolyte  distribution  in  separator 
[11];  the  important  role  of  the  mechanisms  of  oxygen  evolution 
and  recombination  on  thermal  runaway  is  discussed  in  [18-22]; 
thermal  management  systems  have  been  proposed  [4,5,23,24]. 

The  aim  of  the  present  work  was  to  evaluate  the  quasi-steady 
state  characteristics  of  a  VRLA  battery  during  overcharge.  To  achieve 
this  goal  we  measured: 

(a)  the  heat  capacity,  Cp,  and  overall  heat  transfer  coefficient,  ft,  of 
the  battery; 

(b)  the  quasi-stationary  characteristics  of  the  system  and  obtained 
stationary  T  vs.  I  curves  by  performing  a  series  of  stepwise  con¬ 
stant  current  polarization  experiments; 

(c)  the  kinetics  of  the  oxygen  evolution  reaction  proceeding  at  the 
positive  plates  of  VRLA  batteries. 

Based  on  the  obtained  results,  we  discuss  the  conditions  for 
possible  thermal  failure  of  the  battery. 
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2.  Charge  and  discharge  reactions  and  thermal  behavior  of 
VRLA  batteries 

A  VRLA  battery  consists  of  Pb/PbS04  negative  electrodes, 
Pb02/PbS04  positive  electrodes,  absorptive  glass  mat  (AGM)  sep¬ 
arator  and  H2S04  electrolyte.  The  electrochemical  reactions  that 
take  place  at  each  electrode  during  cell  operation  can  be  summa¬ 
rized  as  follows  (left  to  right  are  the  charge  reactions,  and  in  the 


reverse  direction  are  the  discharge  reactions.): 

Positive  electrode  (Pb02/PbS04): 

PbS04(s)  +  2H20  +>  Pb02(s)  +  HS04-  +  3H+  +  2e"  (1 ) 

with  oxygen  evolution 

2H20  -*  02+4H++4e-  (2) 

hydrogen  recombination 

H2  2H+  +  2e_  (3) 

and  grid  corrosion 

Pb  +  2H20  -*  Pb02  +  4H++4e-  (4) 

Negative  electrode  (Pb/PbS04): 

PbS04(s)  +  H+  +  2e"  +>  Pb(s)  +  HSCXp  (5) 

with  oxygen  recombination 

02  +4H+  +4e_  -+  2H20  (6) 

and  hydrogen  evolution 


Table  1 

Thermodynamic  properties  of  VRLA  battery. 


Symbol 

Value 

Unit 

Description 

as2 

+326.412 

JK-1  mol”1 

Entropy  change  of  reaction  (2) 

as6 

-326.412 

J  K-1  mol”1 

Entropy  change  of  reaction  (6) 

as7 

0 

J  K-1  mol”1 

Entropy  change  of  reaction  (7) 

<  (25  °C) 

1.229 

V 

Standard  potential  of  02/H20  at 

25  °C  towards  SHE 

<g/Ag2S04  (25°C) 

0.654 

V 

Standard  potential  of  Ag/Ag2S04  at 
25  °C  towards  SHE 

P/dr)o2/H20 

-0.846 

mVI<_1 

Isothermal  temperature  coefficient 
of  standard  potential  of  02/H20 

( Atp/dl )  Ag/AfcS04 

-1.182 

mVI<_1 

Isothermal  temperature  coefficient 
of  standard  potential  of  Ag/Ag2S04 

During  overcharge  of  the  VRLA  cell,  the  total  heat  generation 
rate,  Qceii.  is  a  sum  of  the  rates  of  the  endothermic  and  exother¬ 
mic  heat  generated  due  to  entropy  changes  of  the  electrochemical 
reactions,  Qs,  the  ohmic  heat  generated  due  to  the  resistance  of 
the  electrolyte  and  the  solid  matrix,  Qohrm  and  the  exothermic  heat 
generated  by  overpotentials  resulting  from  the  shift  from  the  elec¬ 
trochemical  equilibriums  of  the  reactions,  .  The  enthalpy  of  water 
evaporation  is  neglected.  The  sum  +  Qohm  represents  the  Joule 
(irreversible)  heat. 

Qcell  —  Qs  +  Qohm  +  Qjj-  (&) 

As  we  consider  only  the  simplest  electrochemical  mechanisms 
of  02  evolution  and  reduction,  the  above  terms  are  described  as 
follows: 


2H+  +  2e_  H2. 


(7)  (a)  Entropy  term: 


Reactions  (2),  (3),  (6)  and  (7)  constitute  the  internal  oxygen  and 
hydrogen  cycles  in  the  cell.  In  the  presence  of  these  side  reactions, 
the  VRLA  cell  is  a  three-phase  system  consisting  of  the  solid  matrix, 
liquid  electrolyte  and  gas  phase.  At  the  end  of  charge  and  during 
overcharge,  oxygen  is  generated  at  the  Pb02 /electrolyte  interface 
and  may  evolve  in  the  gas  space  when  its  solubility  limit  in  the 
electrolyte  is  exceeded.  The  oxygen  can  then  be  transported,  via 
the  liquid  and  gas  phases,  from  the  positive  to  the  negative  elec¬ 
trode,  where  the  oxygen  gas  may  dissolve  back  in  the  electrolyte 
and  be  reduced  at  the  Pb/electrolyte  interface.  This  process  forms 
the  internal  oxygen  cycle  in  VRLA  batteries. 

Even  though  the  oxygen  evolution  reaction  on  the  positive  plate 
is,  in  general,  a  multi-step  process,  it  can  be  expressed  by  the  overall 
reaction  (2).  The  oxygen  reduction  reaction  on  the  negative  plate, 
however,  is  much  more  complex.  Although  in  most  cases  it  is  faster 
and  predominates  over  the  hydrogen  evolution  reaction,  reactions 
(6)  and  (7)  are  always  in  competition.  A  detailed  description  of  the 
different  possible  mechanisms  of  oxygen  reduction  at  the  negative 
plate  is  given  in  [18-22].  Here,  we  will  consider  only  the  simplest 
mechanism  described  by  Eqs.  (6)  and  (7). 

Hydrogen  is  generated  according  to  reaction  (7)  at  the  negative 
electrode,  when  the  electrode  potential  goes  too  low  (polarization 
becomes  high).  Hydrogen  recombination,  as  described  by  Eq.  (3), 
is  negligible  because  of  its  poor  kinetics  [25,26].  Accumulation  of 
oxygen  and  hydrogen  in  the  gas  phase  contributes  to  pressure  build 
up,  thereby  causing  venting. 

Grid  corrosion  depends  on  electrode  potential  in  a  complex  way. 
If  the  potential  of  the  positive  electrode  rises,  starting  from  its  open 
circuit  value,  grid  corrosion  decreases  in  the  beginning,  reaches  a 
minimum,  and  increases  again  [27-29].  Compared  to  the  OER,  the 
kinetics  of  the  grid  corrosion  reaction  is  much  slower. 

Heat  source  terms :  As  we  are  interested  only  in  the  overcharge 
process,  the  main  charge/discharge  reactions  ( 1 )  and  (5)  do  not  take 
place,  and  hence  only  the  oxygen  and  hydrogen  cycles  are  operative. 


M 


(9) 


which  is  summed  over  the  electrochemical  reactions  proceed¬ 
ing  in  the  cell.  T  is  the  cell  temperature,  AS/  is  the  entropy 
change  of  reaction  i,  F  is  the  Faraday  constant,  n,  is  the  num¬ 
ber  of  electrons  exchanged  by  the  ith  electrochemical  reaction, 
and  7/  stands  for  the  reaction  currents.  In  our  system,  only  three 
electrochemical  reactions  proceed  during  overcharge:  reaction 
(2)  at  the  positive  electrode  and  reactions  (6)  and  (7)  at  the 


negative  one.  Hence, 

h=b2=i 

(10a) 

fe  —  P)2,rec  =  yf 

(10b) 

ii 

X* 

NJ 

II 

1 

X 

(10c) 

y  is  a  measure  for  the  relative  contribution  of  the  oxygen  and 
hydrogen  reactions  to  the  cell  current  at  the  negative  plate. 
It  gives  the  efficiency  of  the  oxygen  cycle.  When  the  oxygen 
recombination  reaction  is  predominant,  y  tends  to  1,  while 
when  hydrogen  evolution  is  operative,  y  decreases.  The  latter 
happens  when  the  negative  plate  becomes  highly  polarized,  in 
which  case  usually  the  cell  is  gassing.  The  values  of  AS,  are  listed 
in  Table  1  [30].  For  calculation  of  AS/  the  value  of  the  entropy 
of  electron,  Se  =  1  /2SH2(g)  =  65.29  J  K-1  mol-1 ,  is  accounted  for 
[31  ].  Substituting  Eqs.  (10a)-(10c)  into  Eq.  (9),  the  entropy  term 
reads: 

Qs  = -T^(l  -  K)J  (11) 

When  the  efficiency  of  the  oxygen  cycle  is  100  %{y=  1 ),  this  term 
is  equal  to  zero. 

(b)  Ohmic  heat  term 

Qohm  —  f  ^el  —  i  —  tye\) 


(12) 
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<p+  and  cp~x  are  the  electrolyte  potentials  at  the  elec¬ 
trolyte/positive  plate  and  electrolyte/negative  plate  interfaces, 
respectively.  This  term  is  always  positive  and  represents  the 
heat  generated  by  the  electrolyte  resistance.  This  term  may 
become  important  at  high  currents  or  low  saturation  levels, 
when  the  resistance  of  inter-electrode  volume  may  increase 
due  to  formation  of  great  number  of  gas  bubbles  inside  the 
cell. 

(c)  Overpotential  term: 

Qji  =  03) 

i 

It  represents  the  heating  due  to  the  electrochemical  reac¬ 
tions  resistances,  rji  are  the  overpotentials  of  the  reactions. 
This  term,  too,  is  always  positive  and  is  often  called  “Joule 
heating  due  to  electrode  resistance”.  The  overpotentials  of  the 
different  reactions  at  the  positive  and  negative  electrodes  are 
expressed  as: 

rj+  =(p+  -(p+-(p°,  =<p~-<p--<p?  (14) 

where  qfi  are  the  equilibrium  potentials  of  different  elec¬ 
trode  reactions  and  <p+  and  (p~  are  the  positive  and  negative 
plate  potentials.  Substituting  the  expressions  for  the  reaction 
currents,  Eq.  (10a)-(10c),  and  for  the  overpotentials  of  the  reac¬ 
tions,  Eq.  (14),  into  Eq.  (13)  the  expression  for  the  overpotential 
heat  source  term  reads: 


Q*=nu-  (?e+l  -  Y)W°o2  -  <P°h2  )]. 


<  =  1.229  V,  cp°2=  OV  (15) 

where  (p^  and  ^  are  the  standard  electrode  potentials  of 
reactions  (2)  and  (3)  against  SHE  and  U=(cp+  -cp~)  is  the  cell 
voltage. 

Substituting  Eqs.  (11 ),  (12)  and  (15)  into  Eq.  (8)  we  obtain  the 
expression  for  the  total  heat  generation  rate: 


Qce  11  —  ^ 


u-o-y)  (Tir+<-<) 


(16) 


Thermal  balance  of  the  cell :  For  a  single  cell  with  immobi¬ 
lized  electrolyte  and  small  thickness  it  can  be  assumed  that  the 
temperature  distribution  inside  the  cell  is  uniform.  Then  the  dif¬ 
ferential  equation  that  describes  the  temperature  changes  can 
be  reduced  to  the  lumped-parameter  model  of  energy  balance, 
e.g.  [32], 


^  (  Cp  T )  +  hAce\\(T  7a)  —  Qcell 


(17) 


where  Cp  is  the  cell  thermal  capacity,  h  is  the  overall  heat  trans¬ 
fer  coefficient,  Acell  is  the  surface  area  of  the  cell  and  Ta  is  the 
ambient  temperature.  The  second  term  on  the  left-hand  side 
represents  the  effect  of  heat  transfer  from  the  battery  to  the 
ambient  air. 

In  this  work  we  will  consider  the  thermal  behavior  of  a  VRLA 
cell  at  steady  state  conditions,  hence  the  time  derivative  term 
in  Eq.  (17)  is  zero  and  the  thermal  balance  is  reduced  to 


frArellC^  Ta)  —  Qcell 


(18) 


Substituting  the  expression  for  the  heat  generation  rate  Eq. 
(16)  into  Eq.  (18)  we  arrive  to  the  following  expression  for  the 
temperature  vs.  current  dependence: 


T-Ta  = 


hAcell 


.  AS2 
~4F~ 


U-O-Y) 


(19) 


Eq.  (19)  shows  that  when  the  efficiency  of  the  oxygen  cycle  is 
100%,  the  energy  that  is  supplied  to  the  system  is  converted  into 


Fig.  1.  Cool-down  curve.  The  characteristic  time  is  extracted  from  an  expo¬ 
nential  fit  and  the  heat  transfer  coefficient  is  calculated:  Mcell  =  0.15JI<-1  s-1 ; 
h  =  5.7JK“1  s”1  m~2. 


heat  and  the  increase  of  the  cell  temperature  depends  only  on 
the  introduced  power  and  the  heat  transfer  coefficient.  At  lower 
efficiencies,  part  of  the  energy  is  removed  with  the  gases  leaving 
the  cell. 


3.  Experimental 

All  experiments  were  performed  with  2  V/4Ah  test  cells,  com¬ 
prising  2  mm  thick  plates  (two  positive  and  three  negative) 
separated  by  3  mm  thick  AGM  separators.  The  dimensions  of  the 
plates  and  separator  were  5.7  cm  x  6  cm,  and  the  dimensions  of 
the  cells  were  6.5  cm  x  3.3  cm  x  8.7  cm.  Separator  sheets  were  not 
placed  between  the  box  walls  and  the  outermost  negative  plates. 
The  acid  concentration  was  1.28  sp.gr.  (5  M).  All  experiments  were 
performed  at  90%  electrolyte  saturation,  at  which  level  the  cell 
operates  in  recombination  regime.  As  the  saturation  level  was  suf¬ 
ficiently  high  and  operation  of  the  closed  oxygen  cycle  in  the  cells 
was  kept  stable,  it  was  hardly  probable  that  dry-out  of  the  outer 
negative  plates  would  occur.  Cell  current,  voltage  and  temperature, 
as  well  as  positive  and  negative  plate  potentials,  were  moni¬ 
tored  during  the  experiments.  The  gas  leaving  the  cells  was  also 
measured. 

3.1.  Cell  heat  capacity  and  heat  transfer  coefficient 

The  heat  capacity  of  the  cell  was  measured  with  the 
help  of  a  self-made  calorimeter  and  the  reported  value  is 
Cp  =  346.7  ±  11.5  JK-1. 

In  order  to  evaluate  the  temperature  changes  it  was  necessary  to 
consider  the  effect  of  heat  transfer  from  the  cell  to  the  surroundings. 
Hence,  the  overall  heat  transfer  coefficient,  h,  had  to  be  measured. 
We  measured  h  for  our  cell  system  directly  by  heating  the  cell 
up  to  70  °C  and  then  cooling  it  by  natural  convection  to  ambient 
temperature.  From  the  cool-down  curve  (Fig.  1)  and  the  value  of 
the  heat  capacity  we  obtained  the  heat  transfer  coefficient  value: 
hAcell  =  0.15  W  K-1 ;  h  =  5.7  W  K-1  nrr2  (the  outer  area  of  our  test  cell 
was  263.16  cm2). 

The  measured  heat  transfer  parameter  is  within  the  range 
reported  in  the  literature  for  free  air  convection.  It  is  assumed  to 
be  constant  throughout  the  experiments.  The  temperature  depen¬ 
dence  on  h  is  neglected  because  the  temperature  difference  does 
not  exceed  40-50  °C. 
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3.2.  Current/ temperature  dependence 

To  measure  the  characteristic  equilibrium  T  vs.  7  curves  of  VRLA 
cells  during  stable  closed  oxygen  cycle  operation  the  following 
experiment  was  performed.  First,  the  cell  was  fully  charged  at  gal- 
vanostatic  regime  /  =  0.2  A  and  its  capacity  was  checked.  Then,  the 
charged  cell  was  put  on  polarization  mode  comprising  step-wise 
current  changes.  Each  polarization  step  continued  until  constant 
values  of  the  monitored  parameters  were  achieved.  The  gassing 
rate,  V0ut*  was  also  monitored  during  the  experiment.  For  the  sake 
of  brevity,  this  method  will  be  referred  to  as  current  pyramid  or 
galvanostatic  pyramid  further  in  the  paper.  The  maximum  applied 
current  was  2  A  =  0.5Ci0. 


3.2 A.  Effect  of  ambient  temperature,  Ta 

Fig.  2  shows  the  changes  in  U,  T,  cp+ ,  cp~  and  gassing  rate  as  a 
function  of  polarization  time  during  galvanostatic  pyramidal  polar¬ 
ization  from  7  =  2  A  to  7 =0.2  A  and  backwards  at  19  °C  ambient 
temperature.  Each  current  change  leads  to  changes  in  all  monitored 
parameters,  which  tend  to  reach  stationary  values.  The  transition 
(non-stationary)  period  for  all  monitored  parameters  lasted  up  to 
3  h. 

The  experimental  data  evidence  that  there  is  a  slight  hysteresis 
in  the  cell  voltage  during  the  downward  and  upward  current  steps. 
But,  in  general,  the  cell  voltage  stays  almost  constant  for  currents 
higher  than  0.4  A.  The  average  cell  voltage  during  the  experiment 
is  77=2.453  V  (Fig.  2a).  The  changes  in  cell  temperature  follow  the 
changes  in  current:  the  temperature  increases  with  increase  off  and 
vice  verse.  The  downward  and  upward  steps  show  that  the  increase 
of  the  cell  temperature  toward  the  ambient  temperature  depends 
only  on  7  (Fig.  2a). 
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Fig.  2.  Galvanostatic  pyramidal  polarization  experiment  at  Ta  =  19°C:  (a)  cell  cur¬ 
rent,  voltage,  and  temperature;  (b)  positive  and  negative  plate  potentials  and  gassing 
rate. 


The  values  of  the  positive  and  negative  plate  potentials  mea¬ 
sured  against  Ag/Ag2S04  reference  electrode  are  plotted  in  Fig.  2b. 
The  positive  plate  potential  depends  on  both  7  and  T.  As  the  exper¬ 
iment  is  carried  out  with  fully  charged  cell,  (p+  is  determined  only 
by  the  kinetics  of  the  OER.  On  the  one  hand,  at  constant  temper¬ 
ature,  the  current  increase  leads  to  increase  of  cp+.  On  the  other 
hand,  at  constant  current,  cp+  decreases  with  increase  of  T.  The 
competition  between  these  two  effects  governs  the  value  of  cp+ 
during  step-wise  current  polarization  (Fig.  2b).  The  value  of  cp+ 
can  be  determined  from  the  Tafel  dependence  for  OER  knowing 
the  kinetic  parameters  (i.e.  activation  energy,  exchange  current  and 
transfer  coefficient).  In  the  following  subsection  we  will  determine 
the  kinetics  of  OER. 

The  polarization  of  the  negative  plate  increases  with  7  (Fig.  2b). 
Flence,  the  overpotential  of  the  hydrogen  evolution  reaction 
increases  and  H2  may  form,  which  may  lead  to  gassing.  The  gassing 
rate  of  the  cell  is  also  plotted  in  Fig.  2b.  It  can  be  seen  from  the 
figure  that,  when  the  absolute  value  of  the  negative  plate  poten¬ 
tial  goes  above  1.05  V  {\(p~\>  1.05  V),  the  cell  starts  gassing.  Flence, 
for  all  currents  higher  than  1A  the  efficiency  of  the  oxygen  cycle  is 
lower  than  100%. 

Fig.  3  presents  the  dependences  of  the  equilibrium  values  of  cell 
temperature  as  a  function  of  cell  current.  The  data  lie  on  a  straight 
line  for  all  steps  for  which  the  cell  is  in  100%  recombination  regime 
(no  gassing). 

The  same  stepwise  experiment  was  performed  at  Ta  =  30  °C  and 
Ta  =  45  °C.  In  order  to  prevent  gassing,  the  maximum  current  was 
chosen  so  as  not  to  exceed  1  A.  The  stationary  values  of  T  vs.  polar¬ 
ization  current  7  are  also  plotted  in  Fig.  3.  The  linear  dependence 
holds  for  all  three  experiments: 

T  =  T0  +  bl  (20) 

where  T0  is  the  temperature  of  the  cell  at  7  =  0.  For  these  experi¬ 
ments  it  corresponds  to  the  ambient  temperature,  T0  =  Ta,  and  b  is 
the  slope  of  the  curve.  The  slopes  of  the  three  straight  lines  are 
almost  the  same:  b  =  15.2KA-1,  16.1  K A-1  and  15.3  KA-1  atl9°C, 
30  °C  and  45  °C,  respectively.  The  very  close  values  of  b  for  the  three 
experiments  indicate  that  the  mechanism  of  heat  generation  within 
the  cell  does  not  change  and  the  temperature  rise  is  a  direct  result 
of  the  current  increase. 

3.2.2.  Effect  of  heat  transfer  to  the  surroundings 

The  experimental  cell  was  put  in  a  thermally  isolated  box 
to  reduce  heat  dissipation  to  the  ambient  air  and  the  same 


Fig.  3.  Dependence  of  cell  temperature  on  applied  current  at  quasi-steady  state 
conditions  in  case  of  free  air  convection. 
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Fig.  4.  Dependence  of  cell  temperature  on  applied  current  at  quasi-steady  state 
conditions:  experiment  in  thermo-box,  Ta  =  25  °C. 


galvanostatic  pyramidal  polarization  experiment  was  performed. 
The  ambient  temperature  was  25  °C  and  the  maximum  applied 
current  was  1  A.  The  heat  transfer  coefficient  of  this  wrapped 
cell  was  experimentally  measured  and  the  obtained  value  was: 
Mceli  =  0.061  WK-1.  The  stationary  temperature/current  curve  is 
plotted  in  Fig.  4  and  the  variations  of  the  monitored  parameters, 
in  Fig.  5.  It  is  evident  from  the  experimental  data  that  the  linear 
temperature/current  dependence  holds  good  for  non-gassing  cur¬ 
rents,  but  the  slope  has  increased  more  than  twice,  b  =  36.2KA-1, 
due  to  restricted  heat  transfer  to  the  ambient  air.  All  monitored 


Time  /  h 


Time  /  h 


Fig.  5.  Galvanostatic  pyramidal  polarization  experiment  of  the  cell  in  thermo-box  at 
Ta  =  25  °C:  (a)  cell  current,  voltage,  and  temperature:  (b)  positive  and  negative  plate 
potentials  and  gassing  rate. 


parameters  exhibit  the  same  behavior  as  for  the  free  air  convection 
case  (previous  subsection). 


3.2.3.  Thermal  characteristics  of  the  cell,  b 

If  we  look  at  the  thermal  balance,  Eq.  (19),  we  can  see  that  the 
coefficient  b  is  equal  to 


u-(i-y)(r^  +  <-<) 


(21) 


b  is  inversely  proportional  to  the  heat  transfer  coefficient.  Obvi¬ 
ously,  if  the  heat  exchange  between  the  cell  and  the  surrounding 
air  is  impaired  ( h  is  lower),  the  value  of  b  will  increase. 
This  is  demonstrated  with  the  experiment  in  the  thermo-box 
(Fig.  4). 

During  a  single  experiment,  the  heat  transfer  coefficient  h  is 
constant.  Flence,  the  linear  T//  dependence  is  obtained  because 
the  group  in  the  brackets  is  constant.  The  entropy  term  in  the 
brackets  depends  on  temperature,  but  for  the  temperature  interval 
of  the  experiment  Te  (20-60  °C),  it  can  be  considered  as  con¬ 
stant:  it  varies  between  0.25  V  and  0.28  V.  In  general  the  cell 
voltage  depends  on  charging  current.  For  our  experiments  (see  Figs. 
2a  and  5a),  the  voltage  changes  with  current  within  a  range  of 
50  mV  and  30  mV,  for  free  air  convection  and  thermo-box  exper¬ 
iments,  respectively.  For  the  two  systems,  hAcell  is  0.15 WK-1 
and  0.061  WK-1,  respectively.  Hence,  a  cell  voltage  change  of 
50  mV  will  lead  to  a  change  of  0.35  °C  in  cell  temperature  in 
the  former  case  and  0.8  °C  in  the  latter.  This  is  within  the 
range  of  the  experimental  error  for  temperature  measurements. 
Hence,  the  temperature  response  of  our  test  cell  is  insensitive  to 
changes  in  cell  voltage  within  50  mV  (even  100  mV)  for  these  high 
polarizations. 

At  100%  efficiency  of  the  oxygen  cycle,  y=  1  and  the  slope 
b  is:  b  =  U/(Mcell).  Then,  using  the  average  cell  voltage  values 
and  the  measured  hAcell  values  for  the  two  types  of  experiments 
(free  air  convection  and  thermo-box),  the  calculated  slopes  are: 
2.453/0.15  =  16.3  K A-1  and  2.36/0.061  =38.7  K A-1,  respectively. 
These  values  are  very  close  to  the  experimentally  obtained  slopes 
(Figs.  3  and  4). 

When  the  cell  is  gassing,  y  =f=  1,  the  second  term  in  the  brack¬ 
ets  (Eq.  (21))  becomes  operative.  Hence,  at  high  currents,  when  U 
remains  constant  and  the  cell  is  gassing  (see  Fig.  2),  the  cell  tem¬ 
perature  is  lower  than  predicted  from  the  linear  dependence  (see 
the  experimental  data  at  Ta  =  19  °C  and  high  currents  in  Fig.  3). 

The  presented  experimental  results  show  that,  when  the  battery 
operates  at  recombination  regime,  the  linear  temperature-current 
dependence  holds  for  all  ambient  temperatures  and  heat  transfer 
coefficients.  Hence,  this  slope  is  a  characteristic  parameter  of  the 
battery  and  it  would  be  more  convenient  to  present  it  in  dimen¬ 
sionless  form: 


5  =  bMcell 
Dnom 


(22) 


where  Hn0m  is  the  nominal  voltage  of  the  battery.  For  the  two 
studied  cases  (free  air  convection  and  thermo-box),  the  value  of 
b  is  ^1.17.  This  value  should  be  approximately  the  same  for  all 
VRLABs.  Hence,  b  can  be  used  for  fast  and  simple  estimation  of 
the  thermal  response  of  the  battery.  One  only  needs  to  know  the 
nominal  voltage  and  /iAcell  of  the  battery.  hAcell  is  easily  measured 
from  the  cool-down  curve  as  shown  above  (see  Fig.  1).  For  exam¬ 
ple,  for  a  12  V  battery  with  Mcell  =  1.5  W  K-1  the  slope  will  be 
b  =  1.17  x  12/1.5  =  9.36  K  A-1 ,  i.e.  a  current  change  of /=  1 A  will  lead 
to  a  temperature  increase  of  9.36  °C. 

The  experimental  data  show  that  the  highest  increase  in  cell 
temperature  is  observed  when  the  cell  is  not  gassing  and  the  current 
is  the  highest  at  a  certain  cell  voltage  U,  i.e.  the  negative  plate  is  com- 
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pletely  depolarized  (< p~  is  close  to  its  open  circuit  value  <p~  =  -1 V 
vs.  Ag/Ag2S04).  At  these  conditions,  the  cell  current  is  determined 
only  by  the  kinetics  of  the  OER  at  the  positive  plate.  The  positive 
plate  overpotential  will  be  approximately  (p+  =  U  -  1.  Hence,  know¬ 
ing  the  value  of  the  parameter  b  for  a  given  cell  and  the  kinetics  of 
the  OER  reaction  we  can  evaluate  the  cell  voltage  as  a  function  of 
ambient  temperature,  which  yields  a  certain  cell  temperature.  To 
discuss  this,  first  we  will  determine  the  kinetic  parameters  of  the 
OER  in  VRLAB. 

3.3.  Kinetics  of  the  oxygen  evolution  reaction  in  recombining 
VRLA  cells 

There  are  numerous  studies  on  the  kinetics  of  the  OER  on  model 
lead  dioxide  electrodes  in  sulfuric  acid  solution  [33-35  and  refs, 
cited  therein].  The  main  question  is  whether  the  kinetic  parame¬ 
ters  of  OER  are  the  same  when  the  reaction  proceeds  on  a  model 
electrode  in  excess  electrolyte  and  on  a  porous  electrode  within  a 
cell  stack.  As  oxygen  evolution  is  a  gas  reaction,  its  parameters  may 
be  affected  by  the  entrapment  of  formed  gas  bubbles  in  the  pores 
of  the  positive  plate  and  separator  when  the  reaction  proceeds  in 
a  VRLA  cell,  which  reduces  the  active  surface  and  contributes  to 
an  uncompensated  resistance  of  the  pores.  For  example,  the  Tafel 
slopes  of  OER  for  plates  reported  in  the  literature  vary  from  75  mV 
to  110  mV  per  decade  [36-38].  Hence,  our  aim  is  to  study  the  kinet¬ 
ics  of  the  OER  that  proceeds  at  the  positive  plate  of  the  test  VRLA 
cell. 

Normally,  only  the  OER  is  expected  to  proceed  at  the  positive 
plate: 

2H20  -*  02  +4H+  +4e- 

The  corrosion  reaction  also  proceeds  therein,  but  its  exchange 
current  is  very  small,  hence,  its  contribution  to  the  total  current 
compared  to  the  current  of  OER  is  considered  negligible. 

At  high  overpotentials  the  reaction  current  Iq2  obeys  Tafel  rela¬ 
tion: 


Fig.  6.  Current/overpotential  curves  at  different  temperatures  for  the  OER  during 
overcharge  of  VRLA  cell. 

slope  at  high  overpotentials.  There  is  no  indicative  dependence  of 
on  temperature  and  it  can  be  assumed  that  aa  02  is  constant. 
The  obtained  average  value  of  aa>o2  is  0.84. 

The  energy  of  activation  and  the  constant  I<  can  be  estimated 
from  the  Arrhenius  plot.  ln(/0)  vs.  T-1.  ln(/0)  is  the  extrapolated 
Tafel  plot  at  cp+  =  cp^ . 

The  standard  potential  of  the  oxygen  reaction  changes  with  tem¬ 
perature,  which  has  to  be  accounted  for  when  the  value  of  ln(/0) 
is  extracted  from  the  current/potential  curves.  Our  reference  elec¬ 
trode  is  Ag/Ag2S04,  whose  standard  potential  against  SHE  also 
changes  with  temperature.  The  simplest  way  to  account  for  this 
dependence  is  [39]: 

^°  =  ^°(25°C)  +  Cr  — 298)^^  (26a) 


Io  2=1  =  Jo  exp 


(23) 


- 

nF 


(26b) 


where  /0  is  the  exchange  current  of  the  OER  at  the  positive  plate, 
is  the  charge  transfer  coefficient,  R  =  8.31  J  mol-1  K-1  is  the  gas 
constant  and  F=  96487  C mol-1  is  Faraday’s  constant.  Eq.  (23)  can 
be  rewritten  in  the  form: 

RT 

r]  =  <P+-<Po ,  =  - - p(lnJ-lnJo),  (24) 

2  aa,02t 

The  exchange  current  depends  on  temperature.  In  the  simplest 
form  this  dependence  can  be  presented  as: 


70  =  I<  exp 


(25) 


where  I<  is  a  coefficient  which  does  not  depend  on  temperature  and 
is  a  characteristic  of  the  system,  and  A E#  is  the  activation  energy 
of  the  process. 

There  are  three  free  parameters  in  the  current/overpotential 
equation:  aa?o2,  K  and  A E#.  To  evaluate  these  parameters  the  test 
cell  was  put  in  a  thermo-stated  water  bath  and  the  I-V  character¬ 
istics  of  the  fully  charged  cell  were  measured  at  15  °C,  20  °C,  25  °C, 
30  °C  and  40  °C.  The  obtained  experimental  data  are  plotted  in  Fig.  6 
(< p+  is  measured  vs.  Ag/Ag2S04  reference  electrode).  It  can  be  seen 
from  the  figure  that,  for  cell  currents  lower  than  0.1  A,  even  though 
the  overpotential  for  OER  is  high,  the  current/potential  curves  start 
to  deviate  from  Tafel  plots.  Table  2  summarizes  the  charge  trans¬ 
fer  coefficients  at  different  temperatures  evaluated  from  the  Tafel 


The  values  of  cp®  and  {d(p9/dT)  for  02/H20  and  Ag/Ag2S04  elec¬ 
trodes  are  also  given  in  Table  1 .  The  temperature  dependence  of  the 
standard  potential  of  OER  against  Ag/Ag2S04  electrode  is  then 

(p%2  =  0.575  +  0.339  x  10"3(T  -  298).  (26a) 

The  Arrhenius  plot  for  the  exchange  current  density  of  oxy¬ 
gen  evolution,  ln(/0)  vs.  T-1,  is  presented  in  Fig.  7.  Full  squares 
give  the  experimental  values  from  the  water  bath  experiment. 
Other  data  points  in  the  graph  are  from  the  pyramidal  polar¬ 
ization  experiments  in  air.  For  the  latter  experiments,  the 
ln(/0)  values  are  calculated  using  the  experimental  values  of 
( p+ ,  T  and  /  for  the  current  steps  at  which  cp+  reaches  equi¬ 
librium,  and  the  measured  transfer  coefficient  aa)02  =  0.84.  As 
expected,  all  data  lie  on  a  uniform  straight  line.  The  values  of 

Table  2 

Parameters  of  the  Tafel  plot  of  the  oxygen  evolution  reaction  at  the  positive  plate  of 
VRLAB. 
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the  other  two  parameters  are  then  extracted  from  this  straight 
line:  ln(JC,  A)  =  34.4  and  E#  =  155.2  kj mol-1.  (The  small  devia¬ 
tions  of  the  experimental  points  for  current  polarization  from 
the  straight  line  for  the  water  bath  experiment  can  be  ascribed 
to  slightly  different  active  areas  of  the  electrodes  in  different 
cells). 
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Fig.  8.  Critical  cell  current  and  voltage  as  a  function  of  Ta  at  T\im  =  60  °C:  (a)  current 
vs.  Ta\  (b)  positive  plate  potential  and  voltage  vs.  Ta. 


The  effective  Tafel  slope  obtained  here  is  between  67  mV  and 
75  mV  per  decade  for  the  temperature  range  15-40  °C. 

4.  Critical  conditions  for  thermal  failure  of  the  battery 

It  is  well  known  that  high  temperatures  are  harmful  for  VRLAB 
operation  because  of  enhanced  corrosion  and  water  evaporation 
rates.  Degradation  of  the  expander  of  the  negative  plate  starts  above 
a  threshold  temperature  of  about  60  °C.  Hence,  we  will  assume  that 
Tlim  =  60oC  is  the  limiting  temperature  for  normal  battery  opera¬ 
tion.  From  the  linear  current  temperature  dependence  (Eq.  (20)) 
and  the  value  of  the  parameter  b  we  can  calculate  the  limiting  cur¬ 
rent  at  which  the  temperature  of  the  cell  reaches  60  °C  as  a  function 
of  ambient  temperature.  The  values  obtained  for  our  test  cell  are 
plotted  in  Fig.  8.  For  a  given  Ta,  the  cell  will  operate  normally  at 
all  currents  below  the  line  in  Fig.  8a,  whereas  at  currents  above 
the  line,  thermal  failure  of  the  battery  will  occur.  At  float  regime, 
the  potential  of  the  positive  plates  of  VRLAB  is  determined  by  the 
kinetics  of  OER,  as  was  discussed  above.  Hence,  for  our  test  cell  we 
can  calculate  the  limiting  positive  plate  potential  at  which  the  cell 
temperature  is  60  °C,  as  a  function  of  Ta,  from  the  limiting  current 
and  the  kinetic  parameters  of  the  OER  reaction.  This  dependence  is 
plotted  in  Fig.  8b.  As  the  worst  case  for  heating  of  the  cell  is  observed 
when  the  negative  plate  potential  is  close  to  the  equilibrium  one 
(open  circuit),  cp+  =  -1 V  vs.  Ag/Ag2S04,  then  the  limiting  safe  polar¬ 
ization  voltage  is  U= cp+  + 1 V.  Critical  cell  voltages  are  also  plotted 
in  Fig.  8b.  At  all  voltages  bellow  the  curve  the  cell  will  operate 
normally,  while  at  voltages  above  the  curve  it  may  suffer  thermal 
failure. 


5.  Conclusions 

The  present  work  shows  that,  at  electrolyte  saturations  of 
around  90%,  the  thermal  behavior  of  the  VRLA  cell  is  determined 
only  by  the  charging  current  and  the  heat  transfer  to  the  sur¬ 
roundings.  For  small  temperature  changes,  in  the  range  of  50-70  °C, 
which  is  the  case  with  VRLA  batteries,  the  overall  heat  transfer 
coefficient  is  constant.  The  experiments  show  that,  for  recombi¬ 
nation  cells,  the  cell  temperature  is  linearly  proportional  to  the 
charging  current.  When  the  cell  is  gassing,  part  of  the  heat  is  dissi¬ 
pated  with  the  gases  that  leave  the  cell  and  the  cell  temperature  is 
lower  than  predicted  from  the  linear  temperature/current  depen¬ 
dence. 

The  dimensionless  slope  b  can  be  used  for  simple  and  fast  esti¬ 
mation  of  the  thermal  response  of  a  VRLA  battery  to  cell  current 
changes. 

On  overcharge,  the  value  of  the  positive  plate  potential  is  deter¬ 
mined  by  the  kinetics  of  the  oxygen  evolution  reaction.  In  the 
present  work,  the  kinetic  parameters  of  OER  proceeding  at  the 
positive  plate  of  VRLA  cell  are  determined.  The  current/potential 
dependence  of  the  OER  obeys  Tafel  relation.  The  experimental 
results  obtained  indicate  that  the  charge  transfer  coefficient  of  the 
OER  proceeding  at  the  positive  plate  of  hermetically  sealed  lead- 
acid  cell  is  constant  and  its  value  is  higher  than  that  measured 
for  Pb02  electrodes  in  H2S04  solution  (the  Tafel  slope  is  smaller 
compared  to  flooded  Pb02  electrodes). 

The  most  critical  conditions  for  self-heating  of  the  VRLA  cell  is 
when  the  negative  plate  is  least  polarized,  i.e.  the  value  of  <p~  is  close 
to  its  equilibrium  (open  circuit)  value.  Otherwise,  the  overpotential 
of  the  hydrogen  evolution  reaction  becomes  high  enough  and  the 
cell  starts  gassing.  Using  the  linear  T//  dependence,  the  cp+  value 
calculated  for  given  I  and  T,  and  the  open  circuit  value  of  cp~,  it 
is  possible  to  predict  the  limiting  charging  float  voltage  for  safe 
operation  of  the  battery. 
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